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Quantum theory is based on a mathematical structure totally different from conventional arith¬ 
metic. Due to the symmetric nature of bosonic particles, annihilation or creation of single particles 
translates a quantum state depending on how many bosons are already in the given quantum system. 
This proportionality results in a variety of non-classical features of quantum mechanics including 
the bosonic commutation relation. The annihilation and creation operations have recently been 
implemented in photonic systems. However, this feature of quantum mechanics does not preclude 
the possibility of realizing conventional arithmetic in quantum systems. We implement conventional 
addition and subtraction of single phonons for a trapped ^^^Yb^ ion in a harmonic potential. In 
order to realize such operations, we apply the transitionless adiabatic passage scheme on the anti- 
Jaynes-Cummings coupling between the internal energy states and external motion states of the 
ion. By performing the operations on superpositions of Fock states, we realize the hybrid compu¬ 
tation of classical arithmetic in quantum parallelism, and show that our operations are useful to 
engineer quantum states. Our single-phonon operations are nearly deterministic and robust against 
parameter changes, enabling handy repetition of the operations independently from the initial state 
of the atomic motion. We demonstrate the transform of a classical state to a nonclassical one of 
highly sub-Poissonian phonon statistics and a Gaussian state to a non-Gaussian state, by applying 
a sequence of the operations. The operations implemented here are the Susskind-Glogower phase 
operators, whose non-commutativity is also demonstrated. 


I. INTRODUCTION 

Due to the symmetric nature of bosons [T] , the proba¬ 
bility of creating or annihilating a boson is proportional 
to the number of bosons in the original state [2]. The 
single boson operations thus bear totally different math¬ 
ematical structure from conventional arithmetic of ad¬ 
dition and subtraction, which causes the predictions of 
quantum mechanics different from those of classical me¬ 
chanics and provides a foundational basis for commuta¬ 
tion relations. Recently, there have been seminal works 
to realize the bosonic operations at the single-boson level 
for the test of foundations and applications of quantum 
mechanics nn. Here, however, we address a funda¬ 
mental question, ’Does the existence of bosonic annihila¬ 
tion and creation exclude the possibility of realizing the 
conventional arithmetic subtraction and addition opera¬ 
tions in quantum mechanics?’ In our work, we answer 
this question by realizing such operations on the phonon 
states of a trapped ion in a harmonic potential. We per¬ 
form these operations on coherent states and superposi¬ 
tions of Fock states, demonstrating a hybrid of conven¬ 
tional arithmetic and quantum parallelism. The seem¬ 
ingly simple operations which match to classical concepts 
of subtraction and addition bring a classical state into a 
non-classical state nearly deterministically. The technol¬ 
ogy we demonstrate here can be used in many quantum 
systems including circuit quantum electrodynamics pH] 
to enable quantum state engineering. 


The conventional addition and subtraction of a particle 
can be written as 

5+ = ^|n + l)(n|,5- = 5]|n-l)(n|. (1) 

n=0 n=l 

where |n) stands for a Fock state of n bosons. takes 
the n-particle state to the (n + l)-state representing an 
addition, while the subtraction operation, S~ , brings the 
n-state to the (n — l)-state. These operations correspond 
to conventional arithmetic which is commonly used in ev¬ 
eryday life but they do not come out naturally in quan¬ 
tum mechanics. Instead, in quantum-mechanics, creation 
and annihilation a operators are introduced, which 
bear the operator relations 

d) = ^ Vn + 1 \n + 1) (n| , d = ^ \/n \n — 1) (n| . 

n=0 n=0 

( 2 ) 

The proportionality factors \/n + 1 and y/n appear due 
to the symmetric indistinguishable nature of bosons [2] . 
The addition or subtraction in quantum domain thus in¬ 
volves the modification of the probability amplitude of 
state due to the excitation n-dependent factor. While the 
quantum operators (§ have been experimentally demon¬ 
strated |3H9], the realization of the conventional opera¬ 
tions 0 is still to be attested in the quantum regime. 
Indeed, the and S~ operators were suggested as the 
elements of a phase operator by Susskind and Glogower 
m- Thus realization of such operations would serve 
as an important stepping stone to study the properties 
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of the Susskind-Glogower phase operator experimentally. 
In this paper, we demonstrate the operations in a near de¬ 
terministic manner, also showing that the technique is a 
resource to create non-Gaussian states efficiently [Ulisl. 
We show that classical states are turned into nonclassical 
ones manifesting highly sub-Poissonian photon statistics 
and negativity in the Wigner function. The versatility of 
the operations for quantum state engineering is demon¬ 
strated by various sequences of the single-phonon opera¬ 
tions. This is contrasted to the bosonic operations real¬ 
ized so far nn. Their success probability is intrinsically 
low (indeed, the higher the fidelity of the operations, the 
lower the success rate); hence a repetition of such the 
operations is practically limited. 

II. EXPERIMENTAL IMPLEMENTATION 

We implement the S~^ and S~ operations on a motional 
mode of frequency ujx for a single trapped ^^^Yb^ ion in 
a three-dimensional harmonic potential [HHiE] as shown 
in Fig. [^a). The harmonic potential is generated by a ro¬ 
tating electric field in the radial axis with trap frequencies 
(cjx, ^y) = (27r)(2.8, 3.2) MHz and a dc voltage in the ax¬ 
ial direction with ujz = (27r)0.6 MHz. The two hyperfine 
states |F = = 0) = |t) and \F = = 0) = \D 

of the Si /2 manifold represent a qubit state with their 
transition frequency cjhf = (27r) 12.6428 GHz. The anti- 
Jaynes-Gummings interaction, i.e., blue sideband tran¬ 
sition, iJaJC = + h.c., is realized by the 

stimulated Raman laser beams with beat-note frequency 
(<^Ri — ^R 2 ) = {^UF + ^x) + A, where Q is the Rabi 
frequency for the qubit transition, r] = Ak^/h/2MuJx 
the Lamb-Dicke parameter, Ak the net wave-vector of 
the Raman laser beams and M the mass of ^^^Yb^ ion. 
The configuration of the laser beams produces the tran¬ 
sition between |i,n) and |t, R + 1) with the oscillation 
frequency of \/n + lr]Cl as shown in Fig. [^b), due to 
the fundamental property of and d operators in Eq. 
§. Therefore, it is fundamentally impossible to transfer 
\i,n) to It, R + 1) in an n independent manner by apply¬ 
ing the blue-sideband beams at a fixed duration of time. 
Fig. 0 c) shows the experimental time-evolution of the 
blue-sideband transition, where the oscillation rates are 
clearly increased by Vn + 1 factor. 

Here we note that an adiabatic transition shown in Fig. 
[Tie), such as the stimulated Raman adiabatic transition 
[17]. will result in removing \/n + 1 and ^/n in the dy¬ 
namics and enable us to realize S~^ and S~ operations. 
We experimentally erase the excitation-dependent factor 
in the population transfer from the |T, n) state to the 
It, R ± 1) state by applying the scheme of transitionless 
driving [TSHIT]. Generally an adiabatic method provides 
a robust transfer from one state to another against seri¬ 
ous fluctuation in control parameters of a system. For a 
perfect adiabatic transfer, however, a relatively slow op- 


(a) system (d) Controls of n. A, p for adiabatic transfer 



(b) Dynamic blue-sideband transition (e) Adiabatic blue-sideband transition 




(c) Blue-sideband Rabi oscillations (f) Experimental Adiabatic transfer 



FIG. 1: Experimental scheme and dynamic and adiabatic 
transition by anti-Jaynes-Cummings interaction, (a) The 
^^^Yb^ system in a harmonic potential. The qubit level in 
Si /2 manifold, \F = 0,?tif = 0) = |^) and \F = l.rup = 0) = 
It) are coupled by the Raman laser beams to produce the 
anti-Jaynes-Cummings interaction or blue-sideband transi¬ 
tion. We denote Q as the Rabi-frequency on the qubit tran¬ 
sition and the A is the frequency difference of the beat-note 
frequency of Raman beams from cc;hf + cjx- (b) The dia¬ 
gram for the standard blue-sideband transition. The Hilbert 
space is composed of the qubit states, |t)and |J.), and phonon 
states of n excitations, |n). The transitions between |J., n) 
and It, n -h 1) would experience evolutions \/n -h 1 times faster 
than the transition between |t, 0) to |t, 1). (c) Probability of 
finding the ion in | t, '^) state as a function of time, which 
clearly manifests + 1-dependence. The arrow at the first 
peak of each oscillation indicates the duration of a 7r-pulse for 
the corresponding transition. The 7r-pulse duration, T^r, of 
the fundamental blue-sideband transition (red) is 13 /xs. The 
dots represent experimental data and solid lines are from the 
fitting to sin^ • (^) adiabatic blue-sideband 

transition whose frequency is independent of motional quan¬ 
tum number n. The Q and A are controlled as the red and 
blue curves. The z/3 in Q is the counter-diabatic term to sup¬ 
press the transition during the evolution. Here Qo = (27r)38.5 
kHz, 13 = 0.075, and Aq = 1.6Qo- (e) The conceptual dia¬ 
gram for the adiabatic blue-sideband transition without the 
+ 1-dependence, (f) The experimental demonstration of 
the adiabatic blue-sideband transitions realized by the tran¬ 
sitionless quantum driving. The total time to execute the 
transitions is 91/xs for any |n), which is about 7 times Tt^. 
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FIG. 2: Schematic diagram and experimental results for the phonon addition, (a) Implementation of addition is composed 
of a TT-pulse of adiabatic blue-sideband transition which transfers all the phonons from |i, n) to |t 5 '^+ 1) uniformly, followed 
by a TT-pulse of carrier transition which couples the internal states resonantly, bringing down to |i, n). This realizes the 

addition operation of Eq. Q. (c) Additions on a superposition state I'lpi) = ^(|n = 0) + |n=l)). The reconstructed 
density matrices, only the real part of them, indicate the hdelity 0.99 (< 0.01) of the initially prepared state and those of 
the final states 0.96(0.01), 0.92(0.01) and 0.87(0.01) after one, two and three times addition, respectively. The purities of the 
output states are 0.92(0.01), 0.81(0.03), and 0.71(0.06), respectively. The numbers in the parentheses represent the sizes of error 
estimated by the maximal-likelihood methods (see supplementary information). This clearly shows the capability of keeping 
coherence, (c) Wigner functions of the n phonon-added coherent state of amplitude a = 0.8 where n = 0,1, 2 and 3. Observed 
negative values in the Wigner function proves the production of non-Gaussian state which results from a simple operation of 
conventional addition. The fidelities are reduced from 0.97(0.01) for the initial state to 0.87(0.01) (one single-phonon addition), 
0.84(0.01)(two additions), 0.85(0.02) (three additions) and purities are changed from 0.99 to 0.93(0.02), 0.93(0.02), 0.80(0.03). 


eration is required, which makes it impractical in many 
quantum systems. We experimentally reduce the dura¬ 
tion of the adiabatic transfer by adding counter-diabatic 
terms that suppress the excitations [ia[2T]. 

In our scheme, similar to the scheme of adiabatic pas¬ 
sage, we sweep the amplitude Q{t) = CIq [sin(7rt/T) + i/3] 
and the detuning A{t) = Aq cos(7rt/T) of the experimen¬ 
tal control parameters in F7aJC: where = (27r)38.5 
kHz, P = 0.075, and Aq = 1.6^2o as shown in Fig. [^d). 
Note that we add the counter-diabatic phase control i/3, 
which makes the state of the system remain in the in¬ 
stantaneous ground state during the fast driving [MU]. 
On top of these controls, we add the compensation term 
for the AC-stark shift due to the coupling of the blue- 
sideband transition to the qubit level (see Appendix A). 
For the cancellation of the dynamic phases differently ac¬ 
quired during the population transfer, we invert the sign 
of Q and reverse the control of A in the middle of the se¬ 
quence. The total duration of T for the transfer is 7 times 
longer than the TT-pulse time of dynamic blue-sideband 
transition on |n = 0). Note that the standard adiabatic 
scheme is about 20 times longer than the TT-pulse dura¬ 
tion for the same fidelity m and a rapid adiabatic pas¬ 


sage needs a coupling strength four times stronger than 
Oo for fho same transfer efficiency [22|. Fig. 0f) clearly 
shows the uniform transfer of population from jj., n) to 
It, n + 1) in the range from n = 0 to 5. 


III. CONVENTIONAL ADDITION OPERATION 

We implement the addition operation in Eq. Q 
by first applying the uniform blue-sideband transfer 
X]n=o It? + 1) (t, n| + h.c. and then TT-pulse of car¬ 
rier transition X]n=o +h.c. as shown in Fig. 

I^a). Our addition scheme deterministically adds one 
phonon independent of the initial phonon number state. 
We observe that quantum coherence is preserved in 
the addition operations. We prepare an initial state 
{\n = 0) + |n = 1)), apply the additions up to three 
times and measure the density matrix of the result¬ 
ing phonon states. As shown in Fig. |^b), the co¬ 
herences represented by the off-diagonal terms of the 
density matrix clearly remain after addition processes. 
We construct the density matrix by using the iterative 
maximum-likelihood algorithm m after displacing the 
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(a) Subtraction operation 

1. Carrier transition 2. Adiabatic blue-sideband transition 




(c) \ii/,) = \a = l2) 



FIG. 3: Schematic diagram and experimental results of 
phonon subtraction. (a) Sequence of subtraction opera¬ 
tions: the sequence of the operations for addition is re¬ 
versed, z.e., a TT-pulse of carrier transition followed by a tt- 
pulse of adiabatic blue-sideband transition, which implements 
l"^’'^) (t?'^ + 1|- This takes the phonon state from \n + 1) 
to \n) except \n — 0) where |^, 0) transfers to |t, 0). The |t, 0) 
state is abandoned by the conditional measurement after the 
fluorescence detection which only collects non-zero phonon 
state data without fluorescence, (b) Subtraction on a su¬ 
perposition state |'0i) = ^ (|n = 2) + |n = 3)). The popula¬ 
tion is reduced and the coherence is conserved. The initial 
hdelity and purity of the state are 0.83(0.02) and 0.73(0.03). 
The fidelities are changed to 0.77(0.02) and 0.83(0.01) after 
one and two times subtraction, respectively. The purities be¬ 
come 0.65(0.02) and 0.75(0.02). In the preparation and dis¬ 
placement operations for the superposition states with flu¬ 
orescent detection, the zero components are increased due 
to unexpected experimental imperfections, which accidentally 
increase the fidelity and purity for the state ^(|0) + |1)). (c) 
Subtraction from an initial coherent state \a — 1.2). The ini¬ 
tial hdelity of the state is 0.96(0.01) and the hdelities are 
reduced to 0.92(0.01) and 0.66(0.01) after one and two times 
subtraction, respectively. 


state by 8 different angles with the amplitude of a ^ 0.8 
and then by hnding the phonon number distributions 
(see Appendix B). The phonon number distributions are 
obtained by observing the time evolutions of the stan¬ 
dard blue-sideband transitions, similar to the direct re¬ 
construction scheme of the phonon density matrix m- 


As a second example, we prepare an initial coherent 
state |a = 0.81) and apply the addition operations as 
shown in Figs. ic) and (d). One immediate conse¬ 
quence of the addition operations on the coherent state 
is the production of sub-Poissonian phonon statistics be¬ 
cause the addition increases the average phonon number 
but not the shape of the distribution and variance. We 
observe that the ratios between the variance and the av¬ 
erage phonon number reduce from 1 to 0.43, 0.39, 0.2, 
after one to three single-phonon additions, respectively. 
Applying the first addition operation, we detect nega¬ 
tivity in the Wigner function as shown in the second 
column of Fig. [^d). It is important to note that the 
addition operation, which converts a coherent state to a 
highly non-Gaussian state, is nearly deterministic unlike 
the case of d) operation m- There is a limit in the 
number of additions we can apply, due to the validity of 
the adiabatic approximation and the heating process of 
phonons [25]. Under this limitation, we could perform 
the operations three times without the significant loss of 
fidelity. As shown in Fig. id), the experimental results 
and the theoretical predictions for the Wigner functions 
are in excellent agreement. This is significant in com¬ 
parison to the photonic realization of bosonic operations 
of single photon creation and annihilation [Tj. Here we 
obtain the Wigner function of the state from the recon¬ 
structed density matrix. 

IV. CONVENTIONAL SUBTRACTION 
OPERATION 

The subtraction operation S~ in Eq. Q is real¬ 
ized by reversing the sequence of the addition oper¬ 
ation, that is, the application of the 7r-pulse of car¬ 
rier transition and the uniform blue-sideband transfer 
X]n=i ~ (t:^l + d.c., followed by the fluorescent 

detection as shown in Fig. [^a). The zero phonon state 
|n = 0) is eliminated after the subtraction, which is im¬ 
plemented by the conditional measurement in our exper¬ 
imental scheme. After the detection sequence, we only 
collect the data with no fluorescence, which has the suc¬ 
cess rate given by the probability of the non-zero phonon 
states. We examine the performance of the subtrac¬ 
tion operation with an initial phonon superposition state 
^ (|n = 2) + |n = 3)). As shown in Fig. I^b), the sub¬ 
traction operation reduces the phonon excitation by one 
quanta. After the second application of the subtraction, 
the off-diagonal terms of the density matrix are signifi¬ 
cantly reduced, which shows the current limit in experi¬ 
ments due to the heating of the system. We also prepare a 
coherence state \a = 1.2) and apply the subtraction twice 
as shown in Figj^c), which shows that qualitatively the 
subtraction works for any initial quantum state. The sub¬ 
traction operation can squeeze a coherent state which is 
different from annihilation that has the coherence state 
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FIG. 4: Experimental results after addition-then-subtraction and subtraction-then-addition, respectively. Only the real part 
of experimentally measured density matrices is shown a. for an initial coherent state |'0i) = \a = 1.2), b. single-phonon added- 
then-subtracted S~I'lpi) state and c. single-phonon subtracted-then-added S~^S~ \^pi) state, b. The state after addition- 
then-subtraction is basically identical to the original state. The fidelity of the S~l^pi) state to the original state I'lpi) is 
0.97(0.01) and the purity is 0.96(0.01). c. The state after subtraction-then-addition is not same as the original state, because 
the vacuum component is thrown away during the projective measurement. The small population in zero component mainly 
comes from the imperfection of the fluorescence detection and heating of the system, which is in good agreement with numerical 
simulation. 


as its eigenstate. However, our experimental precision is 
not high enough to observe the squeezing effect. 

V. COMMUTATION RELATION OF THE 
ADDITION AND THE SUBTRACTION 

We study experimentally how the quantum states are 
changed depending on the order of the addition and sub¬ 
traction for an initial coherent state \^pi) = \a = 1.2). 
If we add then subtract S~ S~^ the state after the 
sequence is the same as the original one, since there is 
no amplitude modification. For the case of subtraction- 
then-addition S~^S~ I'lpi), the final state does not have 
vacuum component because the vacuum state is removed 
at the first subtraction. Fig. [^b) shows the experimen¬ 
tal result of S~S~^ which is basically identical to 
the initial state of Fig. ®a). Fig. |^c) shows the result 
after the operation of S^S~ I'lpi), where there is no sig¬ 
nificant vacuum component in the density matrix. The 
vacuum component is not perfectly removed because of 
the detection error during the projective measurement 
based on the atomic fluorescence and heating of the sys¬ 
tem. The fluorescent detection duration is comparable to 
the motional coherence time of our system, which makes 
the off-diagonal part of the final state suppressed signifi¬ 
cantly (see Appendix C). Our experimental result is well 
in line with the non commuting relation of the Susskind- 
Glogower’s phase operators, i.e., = |0) (0| [TT]. 

VI. REMARKS 

Due to the capability of near-deterministically generat¬ 
ing a non-Gaussian state, the conventional addition and 


subtraction operations provide a new efficient scheme for 
quantum state engineering [26l [SJ. This scheme may 
be further applied to various quantum optics setups such 
as cavity quantum electrodynamics and optomechanics. 
We note a theoretical scheme that the non-Gaussian in¬ 
formation of a phonon state would be transferred to that 
of photon through phonon-photon coupling [28l [29] . It 
has also been discussed that the conventional arithmetic 
addition and subtraction can be used to measure the vac¬ 
uum state without disturbing the state [12] and therefore 
it can be directly applied to construct the Q-function of 
a quantum state of motion. 
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Appendix A: AC Stark shift compensation 

The AC Stark shift in the adiabatic operations mainly 
come from the off-resonant coupling to the carrier tran¬ 
sition, the transition between S 1/2 ^ P 1/2 states of 
^^^Yb^ ion, and the other radial motional mode {uoy ~ 
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ujx + (27r)0.4 MHz). The dominant AC stark shift comes 
from the carrier transition, which is 2 ^^ ^ (27r)33 kHz, 
where Qq = (27r)38.5 kHz and the Lamb-Dicke parame¬ 
ter 77 = 0.89. The amount of the shift brought by Y 
mode is given by that is about 20 times smaller 

than that from the carrier coupling. The AC stark shift 
between qubit states from the Raman laser beams due 

to Si /2 ^ Pi /2 transition is ^^ 2 Ar^^ kHz, 

where gjn and gR 2 are the coupling strengths of Raman 
1 and Raman 2 beam [36] , respectively and the detuning 
from ^Pi/ 2 , = (27r) 18 THz. 


We consider total AC stark shift as the form of , 

and calibrate the blue sideband frequency to be = 

^ ^ pere is the measured blue side¬ 

band frequency. We measure Atotai by varying the in¬ 
tensity of the laser to get several sets of flbsb} 

and by fitting the results with = ^Isb^ + 2 Af^' 

Then we apply the time dependent phase based on the 
amplitude ft{t) = IIq [sin( 7 rt/T) + i/d] and the detuning 
A(t) = Aocos( 7 rt/T) as follows. 


Qo 


sin(7rt/T) cos 


7 ‘ 




— (3 sin 


7 ‘ 


+ A(t'))dt' 


( 3 ) 


Note that the imaginary part here implies | phase differ- based on the measurements of and Atotah 

ence. Here, + A{t'))dt' is calculated as follows 


[\u;l!f+A{t')]dt' = f 
Jo Jo 


, ,meas 
^bsb 


‘^‘AfQfai 

^2 rt 


- Ao cos(7rt/T) 


dt' 


[ [sin^( 7 ri/r) + dt' + AqI sin( 7 rt/T) 

JAfQfal Jo TT 


_ , ,meas4. 

- ^bsb ^ - 




4A 


total 


(1 -h 2/d^)t H-sin(27rt/T) 

27r 


Aq— sin(7rt/T). 

TT 


( 4 ) 

( 5 ) 

( 6 ) 


Appendix B: Reconstruction of phonon density 
matrix 


We use an iterative algorithm proposed in Ref. [23| 
for the reconstruction of an unknown state. It consists of 
a maximum-likelihood estimation solved by expectation- 
maximization algorithm followed by a unitary transfor¬ 
mation of the eigenbasis of the density matrix p. 

Based on the measured phonon distribution of N 
measurements by different displacement, we aim to get 
real probabilities Pn = (n| p|n) that are as close to the 
observed frequencies fn as possible, which can be subject 
to the maximum-likelihood functional 

InL(p) = In JJ(n|p|n)'^” = - '^f„\np„, (7) 

n n 


from which we reconstruct p. This likelihood functional 
can be interpreted as a linear and positive (LP) problem 


in the classical signal processing: 


Pn — 



(8) 


where are eigenvalues of p and hin is a positive ker¬ 
nel. We can solve this LP problem with the expectation- 
maximization algorithm [30l [31] : 


(^) _ ^(^-1) 


= r, 


E 


hinfn 


Pn (r(^-i)) ’ 


(9) 


which is initially set to a positive vector r (r^ > 0 Vi). 

The second part aims at getting the eigenbasis diag¬ 
onalizing the density matrix. This part consists of two 
steps: reconstruction of the eigenvectors of p in a fixed 
basis, and rotation of the basis using a unitary transfor¬ 
mation 

\cf>'j{^’n\ = U\cf>„){K\U^ ( 10 ) 

with the infinitesimal form U = « 1 -h ieG and e 

is a small positive real number. G = i [p, R] is chosen 
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as a Hermitian generator of the unitary transformation, 
where i? is a semipositive definite Hermitian operator 

-R = E„^l«) {n\- 

Starting from some positive initial density matrix p, we 
continue repetition of first finding new eigenvalues Vi us¬ 
ing the expectation-maximization iterative algorithm (§ 
and then finding new eigenvectors (pi by unitarily trans¬ 
forming the old ones. The likelihood of the estimate Pn 
is increased and we finally reach to determine the density 
matrix 


P = '^rn\(pn) {<Pn\- (H) 

n 


Appendix C: Error Analysis 


Dominant error comes from phonon heating process 
caused by the electric-field noise from the trap electrodes 
[32]. Heating decreases the Fock state preparation fi¬ 
delity and affects the adiabatic blue sideband process. 
Its time evolution is known to be described by |32l [33] : 


p{t) = ( 2 d^p {t) a — p (t) da^ — dd^p (t)) + ^ (n + 1 ) ( 2 dp (t) d^ — p {t) dd^ — d^dp (t)) . ( 12 ) 


Here 7 is the coupling strength between the ion motion 
and the thermal reservoir, n is the average phonon of the 
thermal reservoir. In our experimental setup, the heating 
rate yn ~ yn + 1 is 150 Hz. It can be reduced by using 
a large trap, cleaning the electrodes [34| (equivalent to 
reducing y) or cooling the trap (equivalent to reducing 
n). 

More error in experiments may be caused by fitting 
since some noise in the blue-sideband curves may be rec¬ 
ognized as a high phonon population. Although the Fock 
state preparation error will be involved in the transition 
probability, fortunately the error for a small n is rela¬ 
tively insignificant and the population mainly resides in 
small n values for the initial states we prepared in the 
experiments. 

There are some other small errors that we did not cal¬ 
culate or simulate, including the fluctuation of the trap 
frequency ^0.1 kHz coming from ^2% intensity fluctua¬ 
tion of Raman laser, kHz trap frequency fluctuation 
coming from sudden changes of the ion position. We can 
see how much our experimental results are modified when 
comparing to the simulations [35] . 
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